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We study scaling properties of self avoiding polymer stars and networks 
C^- ' of arbitrary given but fixed topology. We use massive field theoretical renor- 

malization group framework to calculate critical exponents governing their 
universal properties (star exponents). Calculations are performed directly in 
J>^, three dimensions, renormalization group functions are obtained in three loop 

approximation. Resulting asymptotic series for star exponents are resummed 
with the help of Pade-Borel and conformal mapping transformation. 
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I. INTRODUCTION 



> 

Many peculiarities in the behavior of polymer chains immersed in a good solvent are 
; well understood qualitatively and described with high precision on the quantitative level 
due to the application of renormalization group methods ]!]-[§ originating from field theory 
B. This progress was initiated by De Gennes ideas linking the model of the <i-dimensional 



Q\ . ferro magnet of 0(n)-symmetry near its critical temperature T = T c to self-avoiding walks 
(SAW) on a lattice of equal dimensionality d |)]||. In particular correlation functions of 
the 0(ra)-symmetric model correspond to those of SVlWs in the limit n = 0. Let us recall 
that the critical exponents v and 7 describing the behavior of the correlation length (£) and 
magnetic susceptibility (x) of the 0(n)-model in the vicinity of critical temperature T c , 



a 

o 

o ■ f ~ r - v 



T — T r 



In polymer theory these exponents correspond to size and configuration number exponents 
v and 7 describing the average square end-to-end distance R 2 

R 2 ^ N 2u , N > 1 (1) 

of a single chain of iV monomers and the number of possible ways of realization of a 
SAW of iV steps on a given lattice: 

Zn^^N^ 1 , iV>l, (2) 

where \i is a non-universal fugacity. 

These laws have been generalized to describe polymers linked at their endpoints to form 
polymer networks |7|-[9|, pd| , p!0|1 . It has been shown that for a given network G (0) holds with 



an exponent 7 = 7g characteristic of the network. Furthermore a subset of networks called 
polymer stars consisting of F polymer chains linked at one endpoint gives rise to a basic 



1 



series of exponents 7f. Just as any network may be disassembled to a number of stars cutting 
polymer chains, its characteristic exponent 7^ is expressed in terms of the exponents 7^ of 
its constituent stars. Whereas the exponents were intensively studied within the 

renormalization group approach to high order and by different methods, this is not the case 
for the star exponents. Their counterparts in the 0(n)-symmetric model were introduced 
as anomalous dimensions of certain composite operators p3j but with no direct physical 
interpretation. In polymer theory they were introduced independently |7|,[5| to characterize 
networks and contacts between chains and later linked to the n = limit of the composite 
operators [I3|. With a number of Monte Carlo simulations at hand P^-jTS| they form 



an attractive and interesting field for analytic calculations. In this article we present our 
results for star exponents obtained in the frames of massive field theory renormalization 
scheme directly in three dimensions. This scheme of calculations as initiated by Parisi fI7 



avoids the £-expansion and thus forms an independent check for the validity of previous 
results [|I^,|T^,|T^,|T^] . As to our knowledge it is applied here in the context of polymer star 
properties for the first time. 

Let us give a short account of the setup of the article. The next section introduces the 
model and defines several series of star exponents giving scaling relations between them. In 
section 3 we calculate the renormalization group functions as series in the coupling to 3rd 
loop order . Being asymptotic these series are resummed in section 4 by means of Pade- 
Borel and conformal mapping technique, we give the numerical values of star exponents 
and discuss them together with results of other approaches. In appendices we show the 
correspondence between the loop integrals entering the theory and their Feynman diagrams, 
and give the details of the pseudo-£-expansion for the star exponents. 



II. POLYMER NETWORKS, STAR POLYMERS AND STAR EXPONENTS 

Consider a polymer network consisting of long polymer chains tied together at their end- 
points (fig. [l]). Mapping such a network on an equivalent field theory one may prove Q |T2|Jl~3 



that as in the case of a linear chain [|T^,|T^] this theory is multiplicatively renormalizable. 
The scaling properties of the network are determined by its "star-like" vertices, connecting 
the extremities of the chains. Let us consider a single star polymer (fig. |^) with F legs of N 
monomers each. It can be shown that for long chains N —>■ 00 the number of self-avoiding 
configurations will scale according to QTUJ : 

2n,f ~ /J FN N^~\ N > 1. (3) 

Formula (j3|) can be considered as the generalization of @, \i still being the non- universal 
connectivity constant and values 7^ give us the first example of star exponents. In the case 
F — 1, 2 one still has the single polymer chain: 71 = 72 = 7 and for F > 3 7^ form a set 
of independent critical exponents. For the two-dimensional case the values of 7^ are known 
exactly 

lF = [68 + 9F(3 - F)]/64. 



This proof is valid for polymer networks of arbitrary architecture at d < 4 and for a large class 
of molecular weight distributions of the strands. 



2 



Exponents 7^ give an example of geometrical exponents associated with special vertices 
where infinite critical objects are fused together. As mentioned in the introduction there is 
no direct analogy to critical exponents describing the 2nd order phase transition in 0(n)- 
symmetric model (see, however formulas (]9j))(^0|) of this article). Let us note here a similar 
problem with an infinite spectrum of scaling dimensions. Consider a d — 2 dimensional 
percolation cluster and the probability Pp of a pinching point joining F peninsulas of the 
infinite incipient cluster. Near the percolation threshold p c this scales like pT|.fZ(1|]: 

Pp~(p-Pcf F , 

where f3p is another basic set of geometrical critical exponents related to 7^ in d = 2. 
Another example are the so-called contact exponents introduced in ]TI| (see [[F] as well) 
to describe the mutual behavior of two polymer stars having F and F' legs as function of 
the distance r of their cores, (fig. [3]). Due to the hard-core repulsion the probability of 
approaching the cores of these two polymers at distance r vanishes as 

P(r)~r 9 ^', r^O. 

Again 9p t pi form a set of universal contact exponents which can not be expressed in terms of 
v and i] ([!]), (H). Let us note that there exist scaling relations between different sets of star 
exponents. Choosing jp exponents as the "basic" ones one gets for the contact exponents 



Of,f> 0: 



vQp^pt — 1 = 7^ + 7f' — Jf+f'- 

Considering a general network G of arbitrary but fixed topology, made of F chains of equal 
length N and tied together in vertices (see fig. [I]) one has for the asymptotic number of 
self- avoiding configurations Zq: 

Z G ~ ^ FN N^°-\ N > 1, 
where the expression for the critical exponent reads: 

7g - 1 + vd = J^np \(j F - 1 + vd) - — (72 - 1 + vd) 

F>1 ^ 

here n F is the number of vertices with F legs. 

Thus knowing one set of exponents for single star polymers one can obtain geometrical 
exponents characterizing any polymer network [@Jl|l2],13|. So in what follows below we will 



consider mainly the properties of a single polymer star immersed in a good solvent with 
arbitrary but fixed number of arms (see fig. |]). We are using the continuous chain model 
introduced by Edwards (see e.g. ||22|| ). Consider an ensemble of F disconnected branches 
and put in correspondence to each branch a a path r a (s), parametrized by < s < S a , a = 
1,2,...,F (S a being the Gaussian surface of ath branch). One can describe such a system 
by the following Hamiltonian: 



p S a 



a=l 

p S a $b 

^ £ JdsJ ds'6 d (r a (s) - r b (s')). 



a, 6=1 
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The partition function is obtained as a functional integral over all possible configurations 
of polymer system divided by its volume Q thus dividing out identical configurations just 
translated in space: 



Z{S a } = yD[r h ...,T F }e W [-^}, (4) 



here the symbol D[ri, rj] includes normalization such that Z{S a } = 1 for w — 0. To 
make the exponential of 5-functions and the functional integral well defined in bare theory 
a cutoff So has to be introduced such that all simultaneous integrals of any variables s and 
s' are cut off by \s — s'\ > so. With this we can relate the Gaussian surface S of a path to 
the notion of steps N of a self avoiding walk by 

N = S/s . 

Introducing into product of 5-functions 

a=2 

ensuring the "star-like" configuration of a set of F chains (c.f. fig. 0), one gets for the 
partition function of the polymer star: 



H 



Z F {S a } = Z F {S 1 ,...,S F } = ^J £>[ri,...,r F ]exp[ 

F 

n^(r a (0) -^(0)). 



F i 

X 



a=2 



Besides the partition functions one can define correlation functions of an F-arm star. 
For the core-endpoint correlation function one has [P,|i^ p~2| , |i~3f : 



P F (r ; r 1; ...,r F] S 1 , S F ) = (f[ 6 d (r a (0) - r )6 d (r a (S a ) - r )}, (5) 



a=l 



here r is the coordinate of the core, ri, being coordinates of chain endpoints, Si, S F 
being their Gaussian surfaces. The averaging in (Bf) means: 



(...) 



i/D[r 1 ,...,i>]exp[-^](. 



Z F {S a } 

For the Green functions one has: 

(27r)^[p + J2 Pa]G F ( Po ; Pi, p F , S x ,..., S F ) 

a=l 



Z F (Si, S F ) / Y[[d d r a exp[ip a r a ]P F (r ;r 1 , ...,r F ; S u S F ). (6) 

a=0 
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The mapping to field theory is performed by a Laplace transformation from the Gaussian 
surfaces S a to conjugated chemical potential variables ("mass variables") ji a : 

f F 

z {Ha} = / II dS b exp[-fi b S b }Z F {S a }. (7) 

J 6=1 

As a result an ensemble of F polymer chains is described by a field theory with the La- 
grangean £{<j) a , /j, a } involving a separate field <p a = a (r) for each chain: 

^M = ^E/^ a |0 a | 2 + ^|V0 a | 2 ) + ^ £ /^r|0 a | 2 |0 a '| 2 , (8) 

a=l a,a'=l 

here 4> a is an n-component vector field (fi a = ((f) 1 , ...,0 n ), \(p a \ 2 = Da=i(0a) 2 - However, the 
perturbation expansion of field theory (|^) at arbitrary value of n results in particular in some 
diagrams, which do not appear in polymer theory . As it is well known [||,|| such diagrams 
(involving closed loops of propagator lines connected to the remainder of the diagram by 
interaction lines only) can be suppressed by taking the limit n — 0. 

Returning now to the partition function of a single polymer star Zp{S b } and using the 
Laplace transformation (0) we can represent its Laplace transform Z F {fi b } in the following 
form: 

ZfM = ^ / d d r n d d r a J D[<p a {r)\ x 

a=l 

F 

II 0a(r o )0a(r a ) exp[-C{(f) b , /x 6 }] | n=0 . 

a=l 

In field theory the star vertex is related to the local composite operator nf=i 4>a of F distinct 
zero component fields. Formally this is the n = limit of an operator known in n-component 
field theory |23|: 

n 

J2 N ai '-' aF ct) ai (r)....(t) aF (r), (9) 

a\ ,...,ap=l 

where N ai, '"' aF is a zero trace symmetric SO(n) tensor: 

n 

jya,a,a s ,...,a F _ g 

a=l 

Now it can be checked diagrammatically that Green functions involving an operator with 
N ai ''"' aF symmetry coupled to F external fields 0™( r ») m ^ ne limit n = coincide with Green 
functions (Q) of the F-arm star polymer. Thus in order to find star exponents in the field- 
theoretical renormalization group program one may consider the problem of renormalization 
of composite operators (Df). This will be the subject of the subsequent section. 
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III. RENORMALIZATION GROUP FUNCTIONS IN THREE LOOP 

APPROXIMATION 



As was shown in the previous section considering the behavior of the star-like vertices 
in polymer theory one faces the problem of renormalization of a field theory containing two 
couplings one being of 0{n) symmetry and described by tensor S ai) ,„ )ai , the other being of 
traceless SO(n) symmetry (jV Q1, " vQF ): 



^N ai '-> aF (f) ai ....(j) aF }. (10) 

Let us point out that the problem of composite field renormalization in the case of traceless 
symmetry we are going to tackle here was considered in the frames of e-expansion in |23f in 
the field-theory context, in polymer context e-expansion results up to e 3 -order were analysed 
in 



[Eq , [L3|l . We will return to these results later. The distinct feature of our study is that 
we are going to apply here equations of the massive field theory. As is well known in order 
to analyze the expressions appearing in such an approach one can either apply e-expansion 
(and thus re-derive the results of ) or as it was proposed by Parisi ]I7| consider 

renormalization group functions directly in the dimension of space of interest (here we are 
interested in the three-dimensional case). So performing calculations in the spirit of [17 
enables us to proceed directly in three dimensions. 

The critical properties of the field theory ( |i~0f) can be extracted from the coefficients 
/3 U , /3 V , J,/,, 7<£2 of the Callan-Symanzik equation for the renormalized iV-point vertex function 

(see e.g. ) . As far as we are interested in the renormalization group functions 

and their derivatives over the coupling constants at the Heisenberg fixed point determined 
by the stable solution of (3 u (u* ^ 0,v* = 0) = we need to calculate vertex functions of 
the symmetric interaction S ai ,...,a 4 , ( j> c ' 1 ■■■■<t> cii F^ 2 \T^ together with F-point vertex function 

with only one N ai, '"' aF (f) c ' 1 ....(f) aF insertion (the other terms will either contain some 
trace of tensor N ai >- ,aF or will be proportional to some power of v and will disappear in 
the Heisenberg fixed point). The graphs for are obtained from the usual four-point 
graphs (see fig. [5]) by replacing each four-point vertex in turn by yN ai, '" ,aF . In three loop 
approximation which we are going to consider here there appear two more graphs in 
which can not be produced in this manner (they are shown in fig. ^). Finally the expressions 
for (3- and 7-functions read 2 : 

(3 U = -(4 -d)u\l-u + (3 2 u LA u 2 + (3 3 u LA u 3 + ..], (11) 



2 We changed the scale for the renormalized couplings u = u (n + 8)D/6, v = v D/6 and beta 
functions (3 u {u) = 6(3 u i (u) /[{n + 8)D], j5 v {u) = 6(3 v '(u)/D, (D being the one-loop integral: D = 
, 2 \ d J jj^iji, u ,v being the renormalized couplings corresponding to the bare couplings uo,vq), 
to define a convenient numerical scale in which the first two coefficients of (3 u {u) are -1 and 1 at 
d = 3. 
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ft = -(4 - d)v[-J^ - — — — — — ^— (n + A 2L V ■+ 



■(4-d) 



n + 8 



0fV + ...], 

2(n + 2) 

l4> = -(4 - ci) . 1 V 2 [2z 2 + (4z 2 - 3* 8 )n + . . .] 
In + o r 



tA J\ n + 2 

7^2 = (4 - d) r w 



n 



l + 7 2LA n + 7 3L V + 



(12) 
(13) 

(14) 



here 8p is the engineering dimension of the coupling v: 

S F 



F + d 

2 



and expressions for two-loop {Pj? LA \ ^ 2LA \ 7 (2LA) ) and three-loop (ffi LA \ (3^ LA \ ^la)) 
contributions read: 



^LA) 

7 (2iA) 
^LA) 

^LA) 



1 
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(3 LA) 



■[(5n + 22)(ii--)+i 2 (n + 2)] 



n + 8) 2LV ~" ~~ /v ~ x 2' 

-(2ii(n + 4F - 2) + 2z 2 (n + 2)/(F - 1) + 

-n-4F + 2))/(n + 8), 
6 



n 



-(1 - 2z 



1 



n + 8) 



i=o 



n 



■[10(n + 8) - (44n + 280)ii + 



+(8 - 3d)(n + 2)z 2 + 12(n + 2)i 3 + 
24(n + 8)i 4 + 6(n + 8)i 5 + 18 (n + 2)i 6 ]. 

In (|15D, (|16D z'o = 1. And for the coefficients one has: 





-8(4n 2 + 61n + 178), 


&i = 


4(31n 2 + 430n + 1240), 




(3dn 2 + 30dn + 48d + 8 




-12(n 2 + 10n + 16), 


bt = 


-48(n 2 + 20n + 60), 


b 5 = 


-24(2n 2 + 21n + 58), 


bl = 


-6(3n 2 + 22n + 56), 


bl = 


-24(5n + 22), 


bl = 


-12(n 2 + 10n + 16); 



(15) 
(16) 



(17) 
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b° v = (n 2 + 8nF + 6n + 20F 2 - 28F + 56), 
bl = -A{n 2 + 7nF + 5n + 18F 2 - 28F + 54), 
b 2 v = -{3dnF - 3dn + 6dF -6d + 4n 2 - 

8nF + 48n - 16F + 80) /(F - 1), 
^ = 12(n + 2), 

b\ = 12(nF + 2n + AF 2 - 10F + 20), 

bl = 3{n 2 + 4nF - 2n + 4F 2 + 12F - 24), 

b 6 v = 6(n + 2F 2 -10F + 18), 

b 7 v = 4(nF - 2n + UF - 28), 

b 8 v = 3(n 2 + 10n + 16)/(F - 1). 

(18) 

here ii — ig are the integrals originating from the corresponding two- and three-loop Feynman 
graphs. Their numerical values at d — 3 are as follows |24|,|25| : 



i x = 2/3; i 2 = -2/27; i 3 = -.0376820725; 
i A = .3835760966; i 5 = .5194312413; i G = 1/2; 

i 7 = .1739006107; i 8 = -.0946514319. (19) 

Correspondence between i\—is and the appropriate Feynman graphs is given in the Appendix 
A. For the additional graphs appearing in one finds the corresponding normalized 
numerical values to be equal: 1 (fig. |]a) and %\ (fig}4|b). 

In the case d = 3 expression for (3 U and gamma-functions coincide with three-loop parts 



of appropriate expressions obtained in ]2B| and at arbitrary value of d they were given in 

Note that expressions (|TTJ) - ([H]) are written for the arbitrary number of field components 
n and thus contain more information that is necessary for the polymer case n = 0. This case 
was considered by Wallace and Zia in three loop ^-expansion motivated by the question 



of relevance of the f-coupling in the renormalization group sense. These authors introduced 
a series of critical exponents otp'- the value of the critical exponent ap is connected with the 
anomalous dimension x' F of v at the 0(n) -symmetrical fixed point: 

F 

(2 - r])a F = x' F - —T], 
, _ dp v (u,v) 

~ dv \u=u*,v=o —Of- l^uj 

The relation to the star exponents is given in the limit n = by 

lF - 1 = - u (2 - rj)a F + [i/(2 - rj) - 1]F. (21) 

Thus the expressions (p|)- ([T^) may be used to study the stability of the 0(n)-symmetric 
fixed point to the perturbation introduced by the traceless coupling. 
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IV. STAR EXPONENTS IN THREE DIMENSIONS 



We now proceed with the calculation of the star exponents. Combining (|Tl"|)-(|T4]) with 
(pip one obtains the following expression for the function 7f(«) (tf = 1f(u*)) 0: 

j F = 1 + (4 - d)F[ 1 1 F LA u + 7 | L V + lf A u% (22) 

lF (23) 

lf A = T—\^K- nF + 7n - + + 10)ii + 
[n + 8)^ 

(dn 2 - dnF + 5dn - 2dF + Qd - 4n 2 + 6nF - 

34n + 8F 2 -4F- 44) /4] , (24) 

7 o = zl(_^ n 3 + rf 2 n 2 F- 
8 

7dV + 4d 2 nF - lQd 2 n + Ad 2 F - lid 2 + 

8dn 3 - Wdn 2 F + 82d * n 2 - 8dnF 2 - 

36dnF + 236dn - 16dF 2 - 32dF + 208d - 16n 3 + 

28n 2 F - 260n 2 + 64nF 2 + 72nF - 1112n + 

80F 3 - 128F 2 + 400F - 1504), 

j F = ^-(dn 2 F - 13dn 2 + 

4dnF 2 + 2dnF - 54dn + 8dF 2 - 56d - 8n 2 F + 100n 2 - 
44nF 2 + 604n - 72F 3 + 152F 2 - 328F + 1000), 

7 | = ^(n 2 -nF + 5n-2F + 6), 

7 3 = -6(-n 2 + nF - 5n + 2F - 6), 

7 £ = _6(-2n 2 + nF 2 + nF - 22n + 4F 3 - 

-14F 2 + 30F-52), 
_3 

7^ = ~Y{n 2 F - 3n 2 + AnF 2 - 6nF - 

18n + 4F 3 + 8F 2 - 36F - 8), 
7 f, = _3(_3 n 2 + nF - I3n + 2F 3 - 

12F 2 + 28F-30), 
1f = -2{nF 2 - 3nF + 2n + UF 2 - 42F + 28). 
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Note scaling relation ([}(]) allowing one to find connection between 7f(^0 and functions (|TT|) 



(14) on the basis of (|20|). Critical exponent n and the combination 2 — v — r] are given by the 
fixed point values of 7^, 7^2 (|l3|) , (|T~^) 
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Calculating star exponents fl22|) one should solve first the fixed point equation: 

(3 u (u*) = (26) 
and then calculate series (|22| ) at u = u* resulting in the value of star exponent 7^: 

7f = 1f ( u = u*). (27) 



In order to be self-consistent substituting result of ( |2T)| ) into the series for 7^ we introduce 
here pseudo-e-expansion (See Appendix B). The resulting series for 7^ ( |B.2|) (as well as 
(0) " (HI)) i s known to be asymptotic and the appropriate resummation procedure is to be 
applied to obtain from it reliable information. 

The results given below were obtained by applying two different resummation techniques 
to the series for 7f(t) ( |B.2| ) as a function of "pseudo-e" parameter r 



7F(r)=Ef/. (28) 
j 

First we tried simple Pade-Borel resummation using [2/1] Pade-approximant 7^(t) for an- 
alytical continuation of the Borel transform 7^(rt) of 7f(t): 

7#(^) = E# (rty 

j 3' 



and writing for the resummed function: 



/ dte-^irt). 
Jo 



Then we applied resummation procedure based on the conformal mapping transformation, 
mapping the domain of analyticity of Jf(t) containing the real positive axis onto a circle 
centered at the origin. Here we introduced a fit parameter (b) considering instead of Borel- 
transform the Borel-Leroy transform Jf~ L (t) of 7f( t ) defined by: 



then 



7 f(t) = / c/tt b e-* 7 |" L (^)- 
Jo 



Assuming that the singularity of 7J? L (r) closest to the origin is located at the point (— 1/a) 
we mapped the r plane onto a circle with a mapping leaving the origin invariant: 

(1 + ar) 1 / 2 - 1 

w = 

' (I + otJVs + i' 



Let us recall that star exponent 7^ equals 7^ = 7f(t = 1). 
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4 w 

T = -- 



a (1 — w) 2 

and thus obtained an expression for ^p~ L {j) convergent in the whole cut plane and, as a 
result, the expression for the resummed function Yf s - While doing this in order to weaken 
a possible singularity on w-plane we multiplied the corresponding expression by (1 — w) a 
and thus introduced one more parameter a. In the resummation procedure the value of a 
(the location of the closest singularity in pseudo-e-expansion) was taken from the known 
large-order behavior of e-expansion series for critical exponents ||28|| , while a was chosen as 



a fit parameter defined by the condition of minimal difference between resummed two- and 
three-loop results. The resummation procedure was quite insensitive to the parameter b 
introduced by the Borel-Leroy transformation (p9|) . 

The results obtained are given in table |j. First we give the values of the exponent 7_p 
obtained by Pade-Borel resummation: the second column contains the value of 7^ obtained 
directly from resummation of the series (|28|), while the third column gives 7^ obtained on 
the base of the resummed series for the exponent x F (namely such a way of calculation of 



7i? was chosen in [|i~3| , |r2f using the e-expansion method). In this case the numerical value of 



'jp was obtained on the base of the numerical value of x F via scaling relation: 

7F = 1 - ux F + [i/(2 - 77/2) - 1]F (30) 

substituting well-known values of the exponents vid = 3) = 0.588, rj(d = 3) = 0.027. 
Columns 4 and 5 give 7^ obtained by the resummation based on conformal mapping tech- 
nique: resummation of the series 7f(t) (fourth column) and resummation of the series x f (t) 
(fifth column). The next columns give the results obtained by the e 3 -expansion based on: 
simple Pade approximation (the 6th column) and Pade-Borel analysis neglecting or exploit- 
ing exact results for d = 2 (the 7th and 8th columns respectively) [13|,12]. The last column 



contains Monte-Carlo data [0,0 ■ For low number of arms F < 5 the results of the different 
approaches agree reasonably well and are also close to the values obtained by MC simulation. 
We have used two different renormalization schemes as well as different procedures for the 
resummation of the resulting asymptotic series. Table [| gives thus a test for the stability of 
the results under changes of the calculational scheme. Obviously for higher number F > 5 of 
arms coincidence of the results is no longer good. The main reason for this is that calculat- 
ing the exponents combinatorial factors lead to an expansion in Fe for the e expansion and 
of Fg when directly expanding in a renormalized coupling g. For such large values of the 
expansion parameter even resummation of the series fails. For large numbers of arms other 
approaches to the theory of polymer stars like a self consistent field approximation might 
be more useful. We conclude that the Parisi method of massive renormalization in fixed 
dimension as it is widely applied in the theory of critical phenomena used together with an 
appropriate resummation scheme is a powerful tool also for the calculation of exponents in 
polymer theory in the present case leading to a good test of previous results and methods. 
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APPENDIX A: APPENDIX 



Here we give the correspondence between loop integrals entering renormalization group 
functions (flTTD - Ql4]) and below) and appropriate Feynman graphs for one-particle irre- 
ducible vertex functions. Fig. |5| give graphs up to the three loop order for the functions T^ 2 ' , 
r^ 4 ) (we keep labeling of ||24j| ). Two additional graphs appearing in r^ F ^ are given by fig]4]a,b. 
In table f| we show the correspondence between numerical values of the loop integrals and 
appropriate Feynman graphs. 



APPENDIX B: APPENDIX 



In the frames of pseudo-e-expansion f\ an auxiliary parameter (let us denote it by r) is 
introduced to keep track of the orders of perturbation theory. The corresponding expression 
for function (3 U ([Til) reads: 



, ) LA nj 

i=2 



P[ LA being j-loop contributions. The equality holds: 

Pu(u,t = 1) = (3 u {u). 
Now one can obtain the value of the fixed point u* as a series in r []: 

u*(r) = Y,u*' {j) T 1 . (B.l) 

3 

This series can be substituted into resulting in the expression for 7f(t) being series in 
powers of "pseudo-5" parameter r. and the final formula for critical exponent 7^ read: 

^ F = Y F eS (u*,T=l). 

Where Yf s m ea n s the resummed (with respect to t) series Jf(t). Performing expansion 
( p.l|) on the base of the expression for the /3-function (|TT|) and substituting the series for the 
Heisenberg fixed point u* 7^ 0, v * = into one finally obtains the following expansion 
of jp in powers of the "pseudo-e" parameter r: 

7*(*) = Et^- (B.2) 

j=0 



5 First introduced by B. Nickel: see Ref.19 in p9|| . 

6 In fact powers of r in a certain term of perturbation theory correspond to the number of loops 
in the loop integrals and this enables one to separate contributions from different orders of the 
perturbation theory while substituting results of ( p6| ) into (p7[). 
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The coefficients 7^ for j = — 3 read: 



7& = 1, 

7^ = (4-rf)F 7 r A) , 

((5n + 22)(i 1 -i) + (n + 2)i 2 ), 
? 3 =(4 _ W+ 16(4-W [(5 „ + 22)( . 1 



n + 8 



(n + 2)z 2 ] + (4 - d)F7j^[^^((5n + 22)(z x - ~) + 

+ {n + 2)i 2 ) 2 + . -— (-32n 2 - 488n - 1424 + 4(31n 2 + 
(n + 8J' 5 

430n + 1240)ii + (3d + 8)(n + 2)(n + 8)z 2 - 
12(n + 2)(n + 8)i 3 - 48(n 2 + 20n + 60)24 - 
24(2n 2 + 2ln + 58)z 5 - 6(3n 2 + 22n + 56)i 6 - 
24(5n + 22)i 7 - 12(n + 2)(n + 8)i 8 )]. 



expressions for 7^ LA , 7j? Lj4 , ^ a are defined in (p3"D - (|2 
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TABLE I. Values of the star exponent jf obtained in three-dimensional theory (columns 2,3,4,5) 
in comparison with the results of e-expansion (columns 6,7,8) and Monte-Carlo simulations (column 
9). See the text for a full description. 



F 



d = 3 
Pade-Borel 



d = 3 
conf.mapping 



e-expansion 



I HP] 



MC 



HIE] 



3 


1.06 


1.05 


1.06 


1.06 


1.05 


1.05 


1.07 


1.09 


4 


0.86 


0.86 


.86 


.83 


0.84 


0.83 


0.85 


0.88 


5 


0.61 


0.61 


.58 


.56 


0.53 


0.52 


0.55 


0.57 


6 


0.32 


0.32 


.24 


.22 


0.14 


0.18 




0.16 


7 


-0.02 


-0.01 


-.17 


-.17 


-0.33 


-0.20 






8 


-0.40 


-0.36 


-.63 


-.62 


-0.88 


-0.60 




(-0.99 


9 


-0.80 


-0.72 


-1.14 


-1.11 


-1.51 


-1.01 







TABLE II. Normalized values of the loop integrals (for graphs 2-M1 - 5-S3 for two-point function 



r( 2 )(fc) the corresponding derivative d/dk 2 |fc2 =0 is given. 



Graph 



Integral 
value 



Graph 



Integral 
value 



Graph 



Integral 
value 



2- U2 

3- U3 

4- U3 

5- U4 

6- U4 

7- U4 



1 
1 

h 
1 

h 
k 



8- U4 

9- U4 

10- U 
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12- U4 
fig.|a 
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k 
1 
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FIG. 1. A polymer network G. It is characterized by the numbers n F of F-leg vertices. Here 
n\ = 3, ri3 = 2, = 1, = 1. 
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FIG. 2. Star polymer. 
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FIG. 3. Two star polymers of functionalities F and F' at a distance r. 
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FIG. 4. Additional graphs appearing in the function T F in three-loop approximation. 
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FIG. 5. Graphs of functions , in three-loop approximation. 
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